
model_yes <- lm(action_yes~condition, data=new_jersey_data)
model_yes_robust <-coeftest(model_yes, vcov.=vcovHC(model_yes, type="HC1")) 
model_yes_rank <- lm(action_yes~condition*rank , data=new_jersey_data)
model_yes_rank_robust <- coeftest(model_yes_rank, vcov.=vcovHC(model_yes_rank , type="HC1")) 
model_yes_trank <- lm(action_yes~condition*t_rank , data=new_jersey_data)
model_yes_trank_robust <- coeftest(model_yes_trank, vcov.=vcovHC(model_yes_trank, type="HC1")) 

stargazer(model_yes, model_yes_rank, model_yes_trank,
		  type = "latex", 
		  digits=2, 
		  star.cutoffs = c(0.05, 0.01, 0.001), 
		  intercept.bottom = FALSE, 
		  se = gaze.coeft(list(model_yes_robust, model_yes_rank_robust, model_yes_trank_robust)), 
		  keep.stat = c("rsq", "n"), 
		  dep.var.labels = "Agreed to Discuss Collaboration",
		  covariate.labels = c("Intercept", "Confidential Condition", "Ranking Condition", "Ranking and Anonymous Condition", "Numeric Rank", "Confidential Condition X Numeric Rank", "Ranking Condition X Numeric Rank", "Ranking and Anonymous Condition X Numeric Rank", "Middle Rank Tercile", "Top Rank Tercile","Confidential Condition X Middle Rank Tercile", "Ranking Condition X Middle Rank Tercile", "Ranking and Anonymous Condition X Middle Rank Tercile","Confidential Condition X Top Rank Tercile", "Ranking Condition X Top Rank Tercile", "Ranking and Anonymous Condition X Top Rank Tercile"), 
		  out = "results/tablee3.tex")
